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Executive  Summary 

Wireless  communications  are  essential  to  emergency  responders.  Unfortunately,  a 
typical  emergency  response  scenario  involves  communication  into  and  within  building 
structures,  which  can  severely  interfere  with  or  completely  block  the  radio-frequency 
channel.  Several  methods  offer  some  potential  relief  to  these  problems,  including  the  use 
of  multi-hop  ad  hoc  networks  for  communication  within  the  building,  and  the  use  of 
smart  antenna  arrays  at  exterior  locations  to  improve  the  communication  between  the 
interior  and  exterior  locations  of  the  building.  However,  due  to  the  enclosed  nature  of 
buildings,  and  the  loss  in  radio-frequency  energy  as  a  signal  propagates  through  floors 
and  walls,  no  single  solution  completely  solves  the  communication  problem. 

In  this  report,  we  consider  another  approach  to  improving  the  communication 
capabilities  of  emergency  responders.  We  investigate  the  possibility  of  creating  antenna 
arrays  from  wireless  devices  such  as  the  emergency  responders'  radios,  wireless  sensor 
networks,  or  arbitrarily  placed  wireless  relay  nodes.  An  antenna  array  allows  the  focusing 
or  directing  of  radiated  electromagnetic  energy  in  a  desired  direction.  By  intelligently 
controlling  the  direction  and  energy  level  of  the  electromagnetic  radiation,  we  can 
increase  the  probability  of  establishing  a  communication  link  and  improve  the  quality  of 
that  link.  The  result  improves  the  emergency  responder's  communication  capabilities. 

Our  approach  is  based  on  antenna  array  theory.  However,  unlike  the  typical 
assumption  that  the  observation  point  or  responder  is  located  in  the  far  field  of  the  array, 
we  assume  the  responder  is  located  within  or  near  the  array  volume.  The  array  volume  is 
essentially  the  building.  We  also  assume  arbitrarily  located  wireless  devices,  which  are  in 
effect  the  antenna  array  elements.  This  initial  investigation  shows  gain  is  possible  from 
arbitrarily  located  wireless  nodes,  at  frequencies  between  1 00  MHz  to  5  GHz.  The  results 
also  indicate  gain  is  possible  whether  the  emergency  responder  is  moving  through  the 
interior  of  the  building,  or  moving  around  the  perimeter  of  the  building.  This  means  the 
array  concept  could  provide  benefit  to  emergency  responders  within  the  building,  and  to 


emergency  responders  communicating  with  people  on  the  exterior  of  the  building.  We 
also  point  out  our  approach  is  not  based  on  ad  hoc  networking  or  protocols.  Although  an 
actual  system  will  make  use  of  such  protocols,  we  are  focusing  on  the  electromagnetic 
behavior  of  the  system  in  aggregate. 

The  preliminary  results  discussed  here  are  encouraging.  We  also  continue  to  improve 
our  model  of  the  real  world  environment  by  refining  the  approximation  used  for  the 
building  materials,  such  as  the  floor  of  the  building.  In  addition,  in  an  actual 
implementation,  emergency  responders  would  likely  quickly  place  wireless  devices  as 
they  enter  the  building,  or  perhaps  even  robots  could  place  the  devices  in  the  near  future. 
We  are  investigating  possible  benefits  from  any  inherent  structure  or  pattern  that  arises 
from  the  placement  of  the  wireless  devices,  or  from  probable  locations  within  the 
building. 
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In  typical  emergency  response  situations  wireless  communications  are  essential  to  the 
emergency  responders.  However,  a  typical  scenario  requires  communicating  into  and 
within  building  structures  that  may  severely  interfere  with  or  completely  block  the  radio- 
frequency  channel.  In  this  report,  we  are  investigating  the  potential  of  utilizing  arbitrarily 
located  wireless  devices  to  increase  the  probability  of  maintaining  a  communication  link 
between  two  points  and  thus  improve  the  emergency  responder's  communication 
capabilities.  The  proposed  approach  employs  methods  based  in  antenna  array  and  general 
optimization  theory.  A  combination  of  analytic  and  simulation  results  allows  rapid  and 
efficient  analysis  of  a  variety  of  array  or  system  configurations.  In  order  to  investigate 
general  trends,  we  approximate  the  floor  of  the  building  as  a  perfect  electrical  conductor 
and  the  wireless  devices  as  Hertzian  dipole  elements.  Observation  points  are  located 
throughout  the  notional  building  volume,  as  well  as  on  a  perimeter  zone.  An  important 
aspect  of  our  analytic  solution  considers  the  observation  points  as  being  in  the  far  field  of 
the  individual  radiating  elements,  with  the  gain  normalized  by  the  system  input  power. 
The  initial  analytic  and  simulation  results  lead  to  the  following  general  conclusions.  First, 
an  appropriately  controlled  system  of  wireless  devices  can  increase  the  communication 
capability  within  a  building  by  optimizing  the  directivity  of  electromagnetic  radiation. 
Second,  the  emergency  responder  can  move  throughout  the  building  and  surrounding 
area,  and  receive  similar  benefit  from  the  array  of  elements.  Third,  the  optimized  current 
phase  information  demonstrates  a  greater  affect  on  the  amount  of  achievable  gain  when 
compared  to  the  optimized  current  magnitude  information. 

Keywords:  Emergency  responder;  random  antenna  array;  optimization;  wireless  devices 
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1.    Introduction 

Wireless  communication  represents  a  key  supporting  technology  to  the  success  of  an 
emergency  responder.  Unfortunately,  a  typical  emergency  response  scenario  involves 
communication  into  building  structures,  which  can  severely  interfere  with  or  completely 
block  the  radio-frequency  channel.  One  potential  method  of  improving  the  radio 
frequency  channel  within  a  building  utilizes  the  intelligent  confrol  of  the  elecfromagnetic 
radiation  from  wireless  devices  quickly  placed  in  the  building  during  entry  by  the 
emergency  responder.  These  devices  would  perform  as  antenna  array  elements  to 
improve  communication  capability  for  emergency  responders  both  within  the  structure 
and  with  external  persons. 

We  here  investigate  the  potential  for  creating  directive  arrays  out  of  arbitrarily  located 
radio  fransceivers.  These  fransceivers  may  take  a  variety  of  forms  such  as  handheld 
radios,  wireless  relay  devices,  or  sensor  nodes.  The  goal  is  to  direct  the  electromagnetic 
energy  for  particular  radio-frequency  channels  within  a  large  volume  area,  (e.g.,  a  large 
building),  so  as  to  improve  both  the  coverage  and/or  data  throughput  of  the  wireless 
communication  channel  between  two  points.  This  research  is  intended  to  directly  improve 
the  design  of  wireless  communication  systems  used  in  current  and  ftiture  emergency 
response  scenarios. 

Our  approach  applies  optimization  techniques  and  antenna  array  theory  to  arbifrarily 
located  radiating  elements.  Optimization  techniques  for  arrays  of  discrete  radiators  have 
a  rich  history,  with  References  [1]  through  [3]  detailing  optimization  techniques  for 
scalar  radiation  from  discrete  radiators  and  chapter  10  in  Reference  [4]  providing  a 
discussion  particularly  relevant  to  this  work.  In  addition,  Reference  [5]  provides  a  highly 
mathematical  approach  to  the  problem  of  optimizing  electromagnetic  radiation  in  a 
particular  direction.  References  [6]  through  [9]  consider  probabilistic  approaches  to 
antenna  array  design,  while  References  [10]  and  [11]  investigate  the  effects  due  to 
random  variation  from  initial  optimized  parameters  such  as  element  location  and  current 


excitation.  However,  to  the  best  of  our  knowledge,  the  appHcation  of  array  theory 
concepts  to  arbitrarily  located  wireless  transceivers  represents  a  new  area  of 
investigation.  In  this  preliminary  report,  we  examine  the  feasibility  of  having  the 
observation  point  at  a  variety  of  locations  within  or  near  the  array.  The  observation  point 
is  analogous  to  an  emergency  responder  trying  to  receive  a  transmission  within  a  building 
that  also  contains  numerous  wireless  devices,  which  constitute  the  array  elements. 
Observation  points  on  the  perimeter  of  the  facility  also  provide  an  important  investigation 
emphasis  since  communication  will  occur  not  only  between  emergency  responders  within 
the  building,  but  with  persons  outside  of  the  building  as  well. 

Our  proposed  topology  of  wireless  devices  resembles  an  ad  hoc  network,  but  our 
focus  is  not  on  ad  hoc  networks  or  protocols.  We  are  studying  the  possibility  of 
controlling  the  electromagnetic  radiation  of  the  wireless  devices  in  aggregate  so  as  to 
form  an  antenna  array.  This  differs  from  the  typical  power  control  of  a  single  device, 
which  is  used  simply  either  to  restrict  or  expand  its  radius  of  coverage,  or  to  increase  the 
battery  life  by  limiting  usage.  A  complete  system  using  our  approach  will  require  various 
supporting  protocols;  however,  our  goal  is  to  determine  the  necessary  electromagnetic 
behavior  of  the  system. 

The  remainder  of  the  report  proceeds  as  follows.  In  Section  2,  we  outline  the  basic 
underlying  theory  behind  our  approach;  in  Section  3,  we  describe  the  general  simulation 
setup  and  scenarios;  in  Section  4,  we  expand  on  further  specifics  for  particular 
simulations  and  draw  important  corresponding  results;  and  finally,  in  Section  5,  we 
present  some  overall  conclusions  and  ideas  for  future  work. 

2.    Basic  Theory 

Basic  concepts  from  antenna  array  theory  serve  as  the  starting  point  for  the  analytical 
basis  of  optimizing  the  electromagnetic  radiation  from  arbitrarily  located  transceivers  or 
elements.  Optimizing  the  ratio  of  two  characteristics  associated  with  the  elements  is  a 


common  approach  used  to  determine  the  appropriate  excitation  and/or  physical 
orientation  of  the  elements.  This  work  seeks  to  maximize  the  electromagnetic  energy  at  a 
given  observation  point  by  determining  the  required  current  (magnitude  and  phase),  on 
each  element.  In  this  work,  we  assume  the  locations  of  the  arbitrarily  located  elements  are 
known.  Figure  1  provides  an  example  of  the  physical  configuration  considered  in  the 
optimization  process. 


Figure  1 .  Example  configuration  with  four  z-directed  dipoles  over  a  ground  plane. 


The  example  configuration  in  figure  1  consists  of  four  z-directed  dipole  elements  over 
a  perfect  electric  conductor  (PEC)  ground  plane.  In  this  analysis,  the  ground  plane  is 
assumed  infinite  in  extent  in  both  the  x  and  y  directions.  The  observation  point  lies  in  the 
vicinity  of  the  elements  ("vicinity"  is  defined  as  lying  within  the  minimal  volume 
containing  the  elements,  or  in  the  perimeter  region  surrounding  the  array).  Two 
characteristics  of  the  configuration  are  used  in  the  optimization  process:  (1)  the  total 
power  provided  to  the  system,  and  (2)  the  z-directed  radiation  intensity  at  the  observation 
point  over  the  ground  plane.  The  optimization  involves  setting  up  the  ratio  of  these  two 


quantities,  and  then  finding  the  complex  current  excitations  that  produce  the  greatest 
ratio. 

Several  typical  measures  of  optimization  apply  to  antenna  arrays,  (e.g.,  see  Reference 
[12]).  However,  we  want  to  optimize  the  electromagnetic  radiation  intensity  for  a 
particular  direction,  i.e.,  the  z-direction,  at  a  specific  observation  point  within  or  near  the 
boundaries  of  the  array  volume.  The  specific  ratio  necessary  to  optimize  takes  the  form 
given  below: 

4^ 

4;r  •  (z  -  component  of  radiation  intensity)  2  rj 

P-  system  input  power  ~  -^-R^^E{r)-T{r)dv 

where  R^^  is  the  average  distance  between  the  elements  and  the  observation  point,  E.  is 
the  z-component  of  the  electric  field,  //is  the  intrinsic  impedance  of  free  space,  E  is  the 

electric  field  vector,  J  is  the  current  vector,  r  is  the  distance  vector  from  the  origin,  and 
dv  is  the  volume  differential.  Distance  i?^^  is  used  because  the  observation  point  is  either 
near  or  within  the  volume  containing  the  array  of  elements,  which  suggests  estimating  the 
distance  to  the  observation  point  by  means  other  than  the  distance  from  the  origin  or 
center  of  the  array  to  the  observation  point.  Note  that  eq  (1)  is  not  exactly  equivalent  to 
the  classical  definition  of  antenna  directivity,  as  given  in  eqs  (2-16)  of  Reference  [13], 
since  that  definition  uses  an  R  value  corresponding  to  the  far  field  of  the  array  and 
assumes  all  the  array  elements  are  located  at  the  origin.  In  our  case,  the  distance  between 
the  observation  points  and  the  elements  is  not  in  the  far  field  of  the  array,  hence 
using /?^^  provides  a  better  mathematical  representation  of  the  overall  behavior  than 
simply  using  some  nominal  R  value,  such  as  the  distance  to  the  farthest  observation  point. 
Our  goal  is  to  maximize  the  value  oip,  which  represents  the  gain  in  a  specific  direction. 

The  array  of  elements  constitutes  a  discrete  system,  and  both  the  numerator  and 
denominator  take  the  form  of  matrices.  In  addition,  the  numerator  matrix  is  Hermitian, 


and  the  denominator  is  Hermitian  positive  definite.  This  allows  the  appHcation  of  classic 
optimization  theory  to  determine  the  maximum  value  ofp  and  the  corresponding  current 
excitations  on  the  elements.  Derivation  of  the  Green's  function,  G%  for  a  z-directed 
Hertzian  dipole  element  over  a  PEC  allows  writing  the  volume  integral  in  closed  form.  In 
particular,  the  electric  field  is  written  as 


where 


G-'iF,FJ='  +;         _     ,  (3) 

'n,  =  ^lix-x„f+{y-y„f  +{z-z„f  ,     (m*  dipole  source  location)  (4) 

m  =  ^|ix  -x„y  +{y-  y^  f  +{z  +  z„  f  ,     (image  of  the  m*  dipole  source)  (5) 


and 


r 


=  V777^,  (6) 


where  F  is  the  distance  from  the  origin  to  the  observation  point.  Other  variables  in  eq  (2) 
are  k,  which  is  the  wavenumber,  j  =  V^ ,  and  a, ,  which  is  the  unit  vector  in  the  z- 
direction.  The  m'  Hertzian  dipole  current  source,  J„,  generating  this  electric  field  is 
given  as 

Jm  =  -a,  — : ,  (/) 

where  co  is  the  radian  frequency,  //q  it  the  permeability  of  free  space,  and  d{) is  the 
impulse  or  delta  function.  By  the  sifting  property  of  the  delta  fiinction  in  the  current 
source,  and  the  dot  product  fi*om  the  vector  multiplication,  the  volume  integral  in  eq  (1) 
takes  the  form 


-\Rc<^E-rdv^UJ-^f^f^G'ir^,FJ  +  ^^'iF^,r„) 


COjU 


0  m=l  n=\ 


k^  dz 


(8) 


Equation  (8)  is  a  square  jVxA'^  matrix,  where  A'^  is  the  number  of  source  elements,  i.e.,  four 
source  elements  means  A'^  =  4.  Let  B  represent  this  denominator  matrix.  Now  letting 
R^  =|F^  -F„|and7?j  =|^m  ~'^n| '  ^^^  taking  the  real  part,  the  terms  in  matrix  eq  (8)  may  be 
written  as 


k        sin(^/?J  ,  sin(^i?,)  ,    1     3^ 


%7tco/iQ       kR^ 


kR, 


k'dz„ 


sm{kRg)     s\n{kR^) 


kR„ 


kR, 


(9) 


Almost  all  the  terms  in  eq  (9)  can  be  written  in  closed  form,  (except  for  the  self  terms 
of  i?^ );  hence,  the  closed  form  expression  for  the  matrix  terms  in  B  is  as  follows: 


B^    = 


v2/  2- 


^  (i+iz„-z„,rk')  ^  3{z„-z„,r 


2d2 


1        3(z„-zJ^ 


2d4 


k'R 


kR„ 


kRi 


v2;,2- 


^     (l  +  (z„+zJ-^0  ,3(z„+zJ- 


2  d2 


1        3(z„+zJ- 
kR,  kRl 


2  n4 


yt'i? 


(10) 


Self  terms  fori?„  =\r^  -^„|  in  eq  (9),  (when  m  =  n),  are  dealt  with  as  a  limiting  case  since 
direct  calculation  leads  to  indeterminate  forms.  Using  a  series  expansion  for  the '— 


terms,  the  self  terms  dependent  on  R^  become 


8;r&jWo 


6  J     ^ttoj/Jq    6 


To  determine  the  closed  form  for  the  numerator  of  eq  (1),  we  take  the  z-component  of 
the  electric  field  derived  in  eq  (2).  Letting R^  =\^ -f^m\  ^"^  ^'m  =\^ ~K\  results  in  the 
following: 


EAr).- 


Ak 


,-Mn. 


kR 


1-- 


kR„ 


{\  +  {z-z„Ye   ^  lAz-zJ"   ^  \z-zj 
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kR'2 


2  D'4 


k'R: 


(11) 


The  m  subscript  on  the  electric  field  refers  to  the  element  or  current  source  responsible 
for  generating  the  field,  as  previously  given  in  eq  (7).  Calling  the  numerator  matrix  A,  we 
can  write  A  as  the  outer  product  of  two  vectors,  i.e.,  A  =  ££* ,  where  e  is  column  vector 
and  f*is  the  complex  conjugate  row  vector.  Using  eq  (11),  the  elements  of  e  are  given 
as 


_     R 


£  = 


.*  Rn 


£„(r)andr=-^£:(F) 


(12) 


We  want  to  maximize  the  ratio    p  =  —  ,  which  equates  to  finding  the  largest 

B 


eigenvalue  of 


Aa  =  AjBa  , 


(13) 


where  a  represents  a  scaling  of  the  current  vector.  Due  to  the  structure  of  the  matrices  A 
and  B,  only  one  nonzero  eigenvalue  exists  for  eq  (13).  We  solve  for  the  maximum  of  this 
ratio  using  the  following  equation 

max{p}=A^^=rB-'6.  (14) 

In  addition,  we  can  determine  the  magnitude  and  phase  of  the  necessary  current 
scaling  fi-om: 

a=B-'£  .  (15) 


Two  important  points  to  note  are  as  follows.  First,  the  scaled  currents  determine  the 
resulting  gain  at  the  observation  point.  Second,  the  existence  of  only  a  single  nonzero 
eigenvalue  follows  from  the  theory  of  eigenvalues  and  eigenvectors.  The  zero 
eigenvalues  correspond  to  eigenvectors  orthogonal  to  the  eigenvector  associated  with  the 
nonzero  eigenvalue.  Physically,  the  currents  corresponding  to  the  zero  eigenvalues 
produce  no  field  at  the  observation  point  [14].  Only  the  currents  associated  with  the 
nonzero  eigenvalue  produce  fields  at  the  observation  point.  A  succinct  development  of  eq 
(13)  through  eq  (15)  is  found  in  Reference  [15]. 

The  sensitivity  of  the  gain  to  the  current  magnitude  and  phase  represents  an  important 

part  of  this  analysis.  Using  the  vector  |a| ,  whose  elements  are  the  magnitudes  of  the 

optimized  currents,  the  sensitivity  of  the  gain  to  the  optimized  current  magnitudes 
follows  from 

«— ;\     \a\  A\a\ 
\a\  B\a\ 

where  D  represents  the  gain  at  the  particular  observation  point  used  in  the  optimization 
process.  In  other  words,  we  simply  calculate  the  gain  using  only  optimized  current 
magnitude  information.  Similarly,  using  Zo ,  the  sensitivity  of  the  gain  to  only  the  phases 
of  the  optimized  currents  follows  from 

D(Ta)=^^  .  (17) 

Za  BZa 

With  the  basic  theory  in  place,  we  now  describe  an  appropriate  simulation  environment  to 
investigate  some  key  attributes  of  the  proposed  approach. 

3.    Simulation  Setup 

Figure  2  below  includes  additional  details  on  the  basic  simulation  configuration 
introduced  in  figure  1.  Three  important  features  in  all  of  the  simulation  runs  include: 


1 .  ground  plane  modeled  as  a  PEC,  infinite  in  the  x  and  y  directions; 

2.  z-directed  Hertzian  dipoles,  with  current  distribution  defined  in  eq  (7); 

3.  observation  points  chosen  in  the  z  =  1 .3  m  plane. 

In  the  simulation  setup,  the  ground  plane  represents  the  floor  of  a  room  or  building.  The 
elements  equate  to  the  wireless  devices  randomly  placed  throughout  the  room  or  building 
as  part  of  the  communication  system  (hereafter  called  simply  the  system).  A  height  of 
1 .3  m  for  the  observation  points  approximates  the  antenna  height  of  the  radio  device  on 
the  emergency  response  personnel,  and  the  z-directed  field  corresponds  to  the  vertically 
polarized  orientation  of  the  antenna.  Figure  2  provides  an  example  of  four  elements  at  a 
height  of  1 .3  m  over  the  ground  plane. 


Figure  2.  Four  z-directed  elements  above  an  infinite  PEC  ground  plane.  All  dipole 
elements  located  at  a  height  of  z  =  1.3  m .  The  footprint  of  the  volume  containing  the 
elements  is  25  m  x  25  m . 


Each  simulation  run  focuses  on  the  effect  of  particular  system  parameters,  which 
include: 

1 .  radio  frequency  range  of  1 00  MHz-5.0  GHz; 

2.  x  ,y  ,  and  z  coordinates  for  observation  points: 
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a.  random  x  and  y  coordinates  uniformly  distributed  over  a  [-8.0m,  8.0m] 
square  (see  figure  3  below); 

b.  deterministic  locations  within  the  array  footprint  (see  figure  4  below); 

c.  random  x  and  y  coordinates  uniformly  distributed  in  a  10  m  perimeter 
around  the  25  m  array  footprint  (see  figure  5  below); 

3.  X  ,y ,  and  z  coordinates  for  element  locations: 

a.  random  x  and  y  coordinates  uniformly  distributed  over  [-12.5  m,  12.5  m] 
square,  height  of  z  =  1 .3  m; 

b.  random  x  and  y  coordinates  uniformly  distributed  over  [-12. 5m,  12. 5m] 
square,  uniform  random  distribution  for  various  ranges  of  z; 

4.  number  of  dipole  elements  ranging  from  2  to  20  elements; 

5.  sensitivity  to  optimized  current  parameters: 

a.  magnitude  sensitivity  calculated  from  eq  (16); 

b.  phase  sensitivity  calculated  from  eq  (17). 
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Figure  3.  Example  of  random  observation  points  within  the  array  in  the  z  =  1.3  m 
plane.  Array  elements  (not  shown)  are  randomly  located  in  either  a  plane  or  volume, 
as  indicated  in  number  3  above. 
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Figure  4.  Deterministic  observation  locations  in  the  z  =  1.3  m  plane.  Array  elements 
(not  shown)  are  randomly  located  in  either  a  plane  or  volume,  as  indicated  in  number  3. 
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Figure  5.  Example  of  random  observation  points  external  to  the  array  elements.  Array 
elements  (not  shown)  are  randomly  located  in  either  a  plane  or  volume,  as  indicated  in 
number  3  above. 

The  simulation  environment  basically  consists  of  using  mathematical  software  to 
randomly  locate  radiating  elements  and  observation  points,  and  then  compute  the  closed 
form  expressions  previously  discussed  in  Section  2.  Each  particular  simulation 
configuration  undergoes  multiple  trials,  e.g.,  1000  trials  with  frequency  =  100  MHz, 
observation  point  x  and  y  coordinates  chosen  at  random,  and  four  dipole  elements. 
Computation  of  the  A  and  B  matrices  follows  from  eq  (10)  and  eq  (12),  while 
optimization  of  the  resulting  gain  is  calculated  from  eq  (14).  Note  in  all  the  simulation 
runs,  the  radiating  or  array  elements  are  located  within  the  notional  building  volume, 
whereas  the  observation  points  may  be  either  interior  or  exterior  to  the  building  volume, 
as  depicted  in  the  three  previous  figures.  Further  details  on  each  set  of  simulations 
precede  the  figures  illustrating  the  corresponding  results  in  Section  4,  along  with  a 
discussion  on  each  set  of  figures.  Section  5  contains  the  overall  conclusions. 
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4.     Simulation  Results 

4.1    Average  Gain  with  Elements  in  a  Plane  and  Interior  Observation  Points 

Each  simulation  run  typically  includes  1000  trials,  (the  only  exception  is  described  in 
Appendix  A).  A  trial  consists  of  first  placing  the  radiating  elements  in  the  notional 
building  space,  then  locating  the  observation  points  in  accordance  with  one  of  the  layouts 
discussed  in  Section  3,  and  finally  calculating  the  gain  for  each  observation  point  by 
optimizing  the  current  magnitude  and  phase  on  each  Hertzian  dipole  element.  Thus,  each 
trial  generates  a  number  of  optimized  gains  equal  to  the  number  of  observation  points, 
(usually  24  in  these  simulations,  with  the  exception  described  in  Appendix  A).  The 
average  gain  is  computed  by  taking  the  average  of  the  average  gain  per  trial  determined 
at  the  observation  points.  A  mathematical  representation  takes  the  form: 

o...=((A))4i:[|t«.l.  ('8) 

where  L  is  the  number  of  trials,  /  is  the  number  of  observation  points  per  trial,  D,  is  the 
gain  at  one  observation  point,  and  Dave  is  the  average  gain.  The  (•)  notation  represents  the 

calculation  of  the  average. 

Figures  6  to  9  depict  simulated  gain  versus  number  of  elements  for  five  different 
frequencies,  with  the  number  of  elements  ranging  from  2  to  20.  Figure  6  shows  the 
average  gain  calculations  based  on  the  x  and  y  coordinates  of  the  observation  points 
uniformly  distributed  over  the  interval  [-8.0m,  8.0m],  (figure  3  depicts  an  example 
distribution).  Figure  7  shows  similar  simulation  results,  but  with  the  deterministic 
observation  points  chosen  to  follow  the  layout  presented  previously  in  figure  4.  Figures  8 
and  9  show  the  standard  deviation  corresponding  to  the  average  gain  of  figures  6  and  7, 
respectively.  This  represents  the  standard  deviation  in  the  average  gain  per  trial  as  given 
by  the  equation 
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a  =  ^{{D„,Y)-{{D^))\  (19) 

where  o  is  the  standard  deviation.  Based  on  this  series  of  plots,  the  method  of  selecting  x 
and  y  observation  coordinates,  such  as  random  versus  periodic,  in  a  specific  plane  of 
observation  does  not  appear  to  significantly  influence  the  statistical  behavior  of  the 
optimization  process. 

Figures  6  through  9  provide  three  general  insights.  First,  the  average  gain 
increases  as  the  number  of  elements  increases.  The  gain  increases  approximately  3  dB 
when  increasing  fi-om  4  to  8  elements  or  from  8  to  16  elements;  or  in  general,  a  doubling 
of  the  elements  causes  a  3  dB  increase  in  gain.  This  is  not  an  obvious  conclusion  since  the 
ratio  under  consideration  as  given  by  eq  (1),  normalizes  the  gain  by  the  total  system  input 
power.  (Note  that  gain  does  not  simply  result  from  an  increasing  number  of  elements; 
compare  to  the  results  in  figures  12,  14,  and  19.)  Second,  the  standard  deviation  appears 
to  settle  below  0.75  dB  in  almost  all  cases.  The  anomaly  for  the  100  MHz  behavior  in 
figures  8  and  9  is  still  under  investigation,  but  we  note  that  the  same  behavior  occurs  for 
both  the  random  and  deterministic  observation  locations. 

The  third  insight  arises  when  we  examine  the  behavior  due  to  the  electrical  height 
of  the  elements  above  the  ground  plane.  Nearly  idenfical  plots  of  the  average  gain  occur 
for  the  frequencies  of  1.0  GHz,  2.5  GHz,  and  5.0  GHz  in  all  four  of  the  figures.  However, 
the  100  MHz  plot  indicates  approximately  an  improvement  of  1  dB  over  these  three 
frequencies  and  a  2  dB  improvement  over  the  500  MHz  case,  for  all  numbers  of 
elements.  This  result  is  due  to  the  height  over  the  PEC  ground  plane  relative  to  a 
wavelength  at  a  specific  frequency  resulting  in  constructive  or  destructive  interference. 
The  use  of  the  PEC  in  the  simulations  represents  a  crude  approximation  to  the  real  world 
surface,  such  as  a  concrete  floor,  so  these  differences  due  to  height  may  be  less 
pronounced  in  the  implementation  of  a  system.  However,  since  the  actual  implementation 
can  be  expected  to  have  devices  located  less  than  2  m  above  the  ground,  the  choice  of 
frequencies  could  significantly  impact  the  actual  gain  achieved. 
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Figure  6.  Average  gain  using  1000  trials  and  24  random  observation  positions  over 
16m  X  16m  square  area. 
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Figure  7.  Average  gain  using  1000  trials  and  24  deterministic  observation  positions  of 
figure  4. 


16 


1.8 

1 1 1 1 1 

- 

1.6 

-e-    100  MHz 

---   500MHz 

m  I'* 

■ 

-+-    1.0  GHz 

- 

2, 

><. 

-B      2.0  GHz 

1  ^-2 

\ 

-V-    5.0  GHz 

- 

CD 

s 

K 

■> 

Q)         1 

■v 

, 

D 

s 

■o 

^ 

1— 

X 

ro   0.8 

[S 

. 

TJ 

/-vV                                                                                                  ^ 

c 

^A 

ro 

^  V  ^                                      "^  , 

^  \    •- 

W    0.6 

"^  0-. 

- 

y^^y__^"---^----o-^- 

^;;-0"--0---0-"^ 

0.4 

~  ^?'"'----s-------_^---_. 

'~  ~*'"~  -  ~  -, 

0.2 
0 

1 1 1 . 1 1 

1              1              1              1              1 

8  10  12  14 

Number  of  Elements 


16 


18 


20 


22 


Figure  8.  Standard  deviation  of  average  gain  plotted  in  figure  6. 
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Figure  9.  Standard  deviation  of  average  gain  plotted  in  figure  7. 
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4.2  Effect  of  Current  Magnitude  on  Optimized  Gain 

Figures  10  through  12  illustrate  the  effects  of  using  only  the  optimized  current 
magnitudes  in  computing  the  gain  results.  Figures  10  and  1 1  result  from  taking  only  the 
magnitude  of  the  optimized  current  vector  for  the  same  configurations  used  to  generate 
figures  6  and  7.  After  determining  the  optimized  current  scaling  from  eq  (15),  application 
of  eq  (16)  provides  the  gain  using  only  the  optimized  current  magnitudes.  All  the 
elements  are  excited  in  phase,  i.e.,  with  a  zero  phase  difference  between  elements. 

As  anticipated  from  the  previous  results,  figures  10  and  11  exhibit  no  significant 
differences.  However,  the  overall  gain  is  substantially  decreased  from  the  gain  depicted 
in  figures  6  and  7.  For  instance,  for  10  elements  in  figures  6  and  7,  all  frequencies  exhibit 
a  gain  in  excess  of  12  dB,  while  for  10  elements  all  frequencies  are  below  8dB  in  both 
figures  1 0  and  1 1 .  Also,  the  rate  of  gain  increase,  approximately  1  dB  per  doubling  of 
elements  in  figures  10  and  11,  represents  a  substantial  decrease  from  the  3dB  per 
doubling  of  elements  depicted  in  figures  6  and  7. 

Figure  12  illustrates  an  important  result  because  it  provides  a  baseline  with  which  to 
compare  current  optimization  techniques.  For  comparison  purposes,  a  random  uniform 
distribution  for  the  current  magnitudes  with  identical  phases  is  used  in  place  of  the 
optimized  current  scaling  vector  in  eq  (16).  The  gain  is  computed  from 

^(^)~.  (20) 

a  Ba 

where  a  is  a  vector  with  a  uniform  random  distribution  for  current  magnitudes  and  zero 
phase.  The  results  are  shown  in  figure  12,  where  the  gain  remains  fairly  constant  as  the 
number  of  elements  increases.  The  gain  in  figures  10  and  1 1  exceeds  the  gain  in  figure  12 
for  all  five  frequencies  and  number  of  elements.  This  indicates  the  optimized  current 
magnitude  information  alone  provides  gain,  even  though  the  amount  of  gain  may  be 
substantially  reduced. 
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Figure  10.  Average  gain  using  only  the  magnitude  of  the  optimized  currents  for  1000  trials 
and  24  random  positions  over  16  m  x  16  m  square  area  in  the  z  =  1 .3  m  plane. 
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Figure  11.  Average  gain  using  only  the  magnitude  of  the  optimized  currents  for  1000  trials 
and  24  deterministic  positions  of  figure  4. 
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Figure  12.  Average  gain  using  uniformly  distributed  random  current  magnitudes  and 
identical  phase  for  1000  trials  and  24  random  positions  over  16  m  x  16  m  square  area  in  the 
z=  1.3  m  plane. 
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4.3  Effect  of  Current  Phase  on  Optimized  Gain 

The  next  two  figures  illustrate  the  effect  of  optimized  phase  on  the  gain.  Figure  13 
illustrates  the  gain  obtained  from  using  only  the  optimized  phase  information  from 
eq  (15)  with  random  observation  points.  This  figure  uses  the  same  points  as  the  average 
gain  plots  in  figures  6  and  10.  These  results  demonstrate  the  significant  effect  of  phase  in 
achieving  the  maximum  possible  gain  for  a  given  configuration.  In  comparing  the  gain 
results  achieved  using  optimized  current  phase  information  only  versus  optimized  current 
magnitude  information  only,  figure  13  versus  figure  10  indicates  a  higher  average  rate  of 
increase  in  gain.  For  example,  at  100  MHz,  the  rate  of  increase  in  gain  is  approximately 
2  dB  per  doubling  of  the  elements  when  using  the  optimized  current  phase  information 
(see  figure  13),  as  compared  to  approximately  1  dB  per  doubling  of  the  elements  when 
using  the  optimized  magnitude  information  (see  figure  1 0). 

Finally,  figure  14  shows  the  results  of  using  a  uniform  random  distribution  for  the 
phase  with  unity  magnitude  for  the  current  vector,  a ,  in  eq  (20).  The  use  of  unity 
magnitude  currents  with  uniformly  distributed  random  phases  over  [0,  27i]  depicted  in 
figure  14  is  nearly  identical  to  the  results  depicted  in  figure  12.  Figures  12  and  14  show 
that  the  optimization  and/or  control  of  the  excitation  currents  (either  current  phase  or 
magnitude),  is  required  to  achieve  gain  at  a  given  observation  point. 
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Figure  13.  Average  gain  using  only  phase  of  optimized  currents  for  1000  trials  and  24 
random  positions  over  16  m  x  16  m  square  in  the  z  =  1.3  m  plane. 
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Figure  14.  Average  gain  using  uniformly  distributed  random  phase  with  unity  magnitude 
currents  for  1 000  trials  and  24  random  positions  over  1 6  m  x  1 6  m  square  area  in  the 
z  =  1.3  m  plane. 
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4.4.  Average  Gain  with  Elements  in  a  Plane  and  Exterior  Observation  Points 

Communication  into  a  building  can  also  be  difficult,  so  the  next  results  focus  on  the 
possible  improvement  in  gain  using  this  multiple  element  system  and  optimization 
process  when  the  observation  point  lies  within  a  10  m  perimeter  of  the  building.  Figure  5 
illustrates  an  example  distribution  of  the  observation  points.  Figures  1 5  and  1 6  depict  the 
average  gain  and  corresponding  standard  deviation  for  the  case  where  the  observation 
point  lies  in  a  10  m  perimeter  of  the  array  boundary.  The  rate  of  gain  increase  in  Figure 
15  indicates  a  gain  of  approximately  3  dB  per  doubling  of  the  number  of  elements,  which 
matches  the  gain  per  element  increase  exhibited  in  figures  6  and  7.  The  general  trend  in 
the  standard  deviations  plotted  in  Figure  1 6  also  mimics  the  trend  shown  in  figures  8  and 
9  (although  o  is  higher  for  some  specific  cases). 

Figures  1 7  and  1 8  show  gain  results  based  on  only  the  optimized  current  magnitudes 
and  phases,  respectively.  Figure  19  depicts  the  use  of  unity  magnitude  current  with 
uniformly  distributed  random  phases  over  [0,  2n]  in  calculating  the  gain.  (The  case  of 
uniform  random  current  magnitudes  with  equal  phases  is  effectively  the  same  as  the 
results  in  figure  19,  and  thus  omitted.)  While  using  either  the  optimized  current 
magnitudes  or  phases  provides  a  gain  increase  over  the  random  current  case,  the 
optimized  phases  enable  a  larger  increase.  The  rate  of  increase  using  the  optimized 
phases  tracks  the  rate  of  approximately  2  dB  per  doubling  of  the  elements  when  using  the 
interior  observation  points.  However,  the  optimized  magnitudes  begin  to  flatten  out  for 
the  14-element  configuration. 
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Figure  15.  Average  gain  for  1000  trials  of  24  random  observation  points  in  the  10m 
perimeter  of  the  array. 
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Figure  16.  Standard  deviation  of  average  gain  for  1000  trials  of  24  random  observation 
points  in  the  10  m  perimeter  of  the  array. 
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Figure  17.  Average  gain  using  only  the  magnitudes  of  the  optimized  currents  for  1000 
trials  and  24  random  observation  points  in  the  10  m  perimeter  of  the  array. 
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Figure  18.  Average  gain  using  only  the  phase  of  the  optimized  currents  for  1000  trials 
and  24  random  observation  points  in  the  10  m  perimeter  of  the  array. 
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Figure  19.  Average  gain  using  uniformly  distributed  random  phase  with  unity  magnitude 
currents  for  1000  trials  and  24  random  observation  points  in  the  10  m  perimeter  of  the 
array. 


4.5  Effect  of  Elements  Located  in  a  Volume  versus  a  Plane 

This  next  series  of  plots  (figures  20  through  31),  depicts  the  gain  behavior  when  the 
radiating  elements  are  not  constrained  to  the  same  plane  as  the  observation  point. 
Specifically,  coordinates  of  the  radiating  elements  are  uniformly  distributed  along  the 
three  coordinate  axes,  x,  y,  and  z,  while  the  observation  points  all  lie  in  the  z  =  1.3  m 
plane.  The  first  six  plots  illustrate  the  gain  and  standard  deviation  for  three  different 
volumes  with  observation  points  within  the  building  boundaries,  while  the  next  six  plots 
result  from  observation  points  in  the  10  m  perimeter.  In  all  six  cases  (three  each  for  the 
interior  and  exterior  observation  points),  the  x  and  y  coordinates  of  the  elements  are 
random  uniformly  distributed  in  the  interval  [-12.5  m,  12.5  m].  The  z  coordinates  are 
random   uniformly   distributed   as   well,   over   one   of  the   three   possible   intervals. 
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[0.1m, 0.5m],  [0.1m,  1.0m],  or  [0.1m,  1.5m].  Note,  these  volume  intervals  start  at 
z  =  0. 1  m  in  anticipation  of  the  elements  being  some  finite  height  above  the  ground  in  a 
practical  system. 

In  the  first  three  plots  (figures  20,  21,  and  22),  the  average  gain  does  not  change 
significantly  as  the  volume  of  possible  element  locations  is  increased.  Also,  the  gain 
exhibits  the  same  rate  of  increase  as  when  the  elements  all  lie  in  the  z  =  1.3  m  plane, 
giving  approximately  3  dB  per  doubling  of  the  elements  (compare  to  figures  6  and  7). 
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Figure  20.  Average  gain  using  1 000  trials  and  24  random  observation  positions  over 

16  m  X  16  m  square  area;  random  radiating  elements  in  a  volume 

25  m   X  25  m   x  0.4  m;  z  positions  range  over  the  interval  [0. 1  m,  0.5  m]. 
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Figure  21.  Average  gain  using  1000  trials  and  24  random  observation  positions  over 

16m  X  1 6  m  square  area;  random  radiating  elements  in  a  volume 

25  m  X  25  m   x  0.9  m;  z  positions  range  over  the  interval  [0. 1  m,  1 .0  m]. 
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Figure  22.  Average  gain  using  1000  trials  and  24  random  observation  positions  over 

16  m  X  16  m  square  area;  random  radiation  elements  in  a  volume 

25  m   X  25  m   x    1 .4  m;  z  positions  range  over  the  interval  [0. 1  m,  1.5  m]. 


These  next  three  plots  (figures  23,  24,  and  25),  illustrate  the  standard  deviation 
associated  with  the  three  previous  gain  plots.  In  all  three  case,  the  standard  deviation  falls 
below  1  dB  for  all  frequencies  after  only  six  elements.  Placing  the  elements  in  a  volume, 
as  opposed  to  a  single  plane  appears  to  smooth  the  standard  deviation.  The  combination 
of  the  gain  and  standard  deviation  indicate  an  overall  smoothing  effect  from  allowing  the 
elements  to  reside  at  a  variety  of  heights.  An  actual  system  will  likely  utilize  elements  at 
a  variety  of  heights,  so  this  smoothing  behavior  without  sacrificing  gain  represents  an 
added  benefit. 


29 


1.8 

1 — 

+ 

- 

-©-    100  MHz 

1.6 

---    500MHz 

- 

' 

-+-    1.0  GHz 

m  1-4 

\ 

^ 

-  B-   2.0  GHz 

. 

■a 

>  ^ 

-V-   5.0  GHz 

c   .  _ 

v\ 

1  ^ 

O    1.2 

\      ^ 

- 

*4_r 

ni 

^ 

^'      ^v 

> 

"*            ^ 

S      ^ 

■ 

^             \><^    ^^ 

- 

-g 

\    °.       -.. 

CD     0.8 

. 

\        N        + :: . 

. 

-o 

^                             -^                              ■-    V 

c 

1      "■□.      ^  =  ,^ 

W    0.6 

0.4 

0.2 

n 

o 

1            1            1           1           1           1           1            1            1 

- 

8  10  12  14 

Number  of  Elements 


16 


18 


20 


22 


Figure  23.  Standard  deviation  of  average  gain  results  shown  in  figure  20. 
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Figure  24.  Standard  deviation  of  average  gain  results  shown  in  figure  21. 
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Figure  25.  Standard  deviation  of  average  gain  results  shown  in  figure  22. 


The  next  series  of  plots  (figures  26,  27,  and  28),  shows  the  average  gain  and  the 
corresponding  standard  deviation  when  the  observation  point  resides  in  the  10  m 
perimeter  of  the  element  volume.  All  the  observation  points  lie  in  the  z  =  1.3  m  plane  and 
are  randomly  distributed.  The  behavior  is  very  similar  to  the  case  with  interior 
observation  points,  namely  a  smoothing  effect  on  the  standard  deviation  while  still 
maintaining  the  same  rate  of  gain  increase  as  when  the  elements  all  reside  in  the  z  =  1 .3  m 
plane.  The  spread  of  the  average  gain  between  frequencies  decreases  from  approximately 
6dB  in  the  planar  case  to  less  than  4dB  in  the  volume  case.  As  mentioned  for  the 
previous  six  plots,  this  indicates  a  potential  added  benefit  for  an  actual  system. 
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Figure  26.  Average  gain  using  1000  trials  and  24  random  observation  positions 

over  10  m  perimeter;  random  radiating  elements  in  a  volume 

25  m   X   25  m   x  0.4  m;  z  positions  range  over  the  interval  [0.1  m,  0.5  m]. 
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Figure  27.  Average  gain  using  1000  trials  and  24  random  observation  positions 

over  1 0  m  perimeter;  random  radiating  elements  in  a  volume 

25  m   X  25  m   x  0.9  m;  z  positions  range  over  the  interval  [0.1  m,  1.0  m]. 
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Figure  28.  Average  gain  using  1000  trials  and  24  random  observation  positions 

over  1 0   m  perimeter  random  radiation  elements  in  a  volume 

25  m   X  25  m   x    1.4  m;  z  positions  range  over  the  interval  [0.1  m,  1.5  m]. 


These  next  three  plots  (figures  29,  30,  and  31),  illustrate  the  standard  deviation 
associated  with  the  three  previous  gain  plots.  In  all  three  plots,  the  standard  deviation 
drops  below  1  dB  for  greater  than  four  elements.  The  spread  between  frequencies  is 
slightly  greater  than  when  the  observation  points  reside  on  the  interior  of  the  building,  but 
the  difference  is  not  large. 
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Figure  29.  Standard  deviation  of  average  gain  results  shown  in  figure  26. 
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Figure  30.  Standard  deviation  of  average  gain  results  shown  in  figure  27. 
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Figure  31.  Standard  deviation  of  average  gain  results  shown  in  figure  28. 
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5.    Conclusions  and  Future  Work 

Several  preliminary  conclusions  related  to  communication  systems  for  emergency 
responders  may  be  drawn  from  these  initial  analytic  and  simulation  results. 

>  An  appropriately  controlled  system  of  wireless  devices  can  theoretically 
increase  the  communication  capability  within  a  building  by  effectively 
functioning  as  an  antenna  array  to  support  optimized  directivity  of  the 
electromagnetic  radiation.  This  can  improve  the  probability  of  receiving  a 
radio-frequency  signal,  as  well  as  the  quality  of  the  signal  in  terms  of  signal- 
to-noise  ratio. 

>  The  emergency  responder  can  move  throughout  the  building  and  surrounding 
area,  and  receive  a  similar  benefit  from  the  array. 

>  The  average  gain  increases  by  approximately  3  dB  per  doubling  of  the 
elements  when  using  optimized  current  magnitudes  and  phases.  In  addition, 
the  overall  average  gain  is  greater  than  14  dB,  with  a  standard  deviation  of  less 
than  0.75  dB  for  all  cases  involving  20  elements. 

>  The  use  of  the  optimized  current  phase  information  only  provides 
approximately  2  dB  of  gain  increase  per  doubling  of  the  elements.  On  the 
other  hand,  using  only  the  optimized  magnitudes  provides  1  dB  or  less  of  gain 
per  doubling  of  the  number  of  radiating  elements.  This  implies  the  achievable 
gain  is  affected  more  by  the  accuracy  of  the  optimized  phases  than  the 
optimized  magnitudes. 

>  Placing  the  elements  in  a  volume  versus  a  single  plane  may  actually  benefit  a 
system  utilizing  several  frequencies  by  smoothing  out  the  gain  statistics. 

One  of  the  challenges  of  implementing  a  system  based  on  this  optimization  process  is 
the  control  of  the  relative  phase  between  elements.  An  important  ongoing  investigation 
examines  the  effect  on  gain  as  the  phases  deviate  from  the  optimal  case.  As  noted  in  the 
fourth  bullet  above,  gain  is  still  achieved  when  only  the  optimized  current  magnitudes  are 
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used,  so  even  if  the  phase  cannot  be  tightly  controlled,  gain  should  be  achievable  since 
the  current  magnitude  is  fairly  straightforward  to  control. 

These  preliminary  results  are  encouraging,  but  further  refinement  of  the  simulation 
environment  is  necessary.  One  improvement  currently  under  investigation  is  the 
modeling  of  the  ground  plane  as  other  than  a  PEC.  Another  point  of  investigation  is  the 
effect  of  mutual  coupling  between  radiating  elements.  Location  sensitivity  analysis  will 
help  determine  how  often  the  optimization  process  would  require  updating  in  reaction  to 
changes  in  the  topology  of  the  system  of  wireless  devices.  Finally,  the  impact  of  losing 
some  of  the  array  elements  must  be  studied,  as  the  actual  practical  implementation  will 
operate  in  a  harsh  environment. 
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Appendix  A:  Examination  of  Statistics  Gathering  Metliod 

This  series  of  plots  demonstrates  the  effect  of  choosing  1 000  observation  points  and 
24  trials  versus  24  observation  points  and  1000  trials.  Figures  A-1  and  A-2  depict  the 
average  optimized  gain  (i.e.,  using  optimized  magnitude  and  phase  for  element  currents 
to  achieve  gain  at  a  specific  observation  point,  as  performed  in  Section  4.1),  when  the 
observation  points  reside  in  the  interior  and  the  exterior  of  the  building,  respectively.  A 
uniform  random  distribution  on  both  the  x-  and  y-axis,  and  the  z  =  1.3  m  plane  provide 
the  observation  point  coordinates.  Compare  these  first  two  plots  to  figures  7  and  15, 
respectively.  The  results  are  nearly  identical  between  figure  A-1  and  figure  7,  and  figure 
1 5  and  figure  A-2. 

Figures  A-3  and  A-4  illustrate  the  corresponding  standard  deviation  of  the  average 
gain  in  figures  A-1  and  A-2.  These  two  plots  compare  to  figures  9  and  16.  All  the  trends 
are  similar,  including  the  fi^equency,  which  shows  the  greatest  standard  deviation,  namely 
500  MHz  for  figures  34  and  9,  and  1.0  GHz  for  figures  A-3  and  16.  Since  computing  the 
standard  deviation  of  the  average  per  trial  generates  these  plots,  some  differences  are 
expected.  However,  statistically,  the  results  are  almost  identical,  which  implies  either 
method  of  collecting  first-order  statistics  will  provide  similar  results. 
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Figure  32.  Average  gain  using  24  trials  and  1000  random  observation  positions 
over  16  m  X  16  m  square  area  in  the  z  =  1.3  m  plane. 
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Figure  33.  Average  gain  for  24  trials  of  1000  random  observation  points  in  the 
10  m  perimeter  of  the  array. 
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Figure  34.  Standard  deviation  for  average  gain  results  in  figure  32. 
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Figure  35.  Standard  deviation  of  average  gain  results  shown  in  figure  33. 
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